Quantum charge glasses of itinerant fermions with cavity-mediated long-range interactions 
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We study models of itinerant spinless fermions with random long-range interactions. We motivate such 
models from descriptions of fermionic atoms in multi-mode optical cavities. The solution of an infinite-range 
model yields a metallic phase which has glassy charge dynamics, and a localized glass phase with suppressed 
density of states at low energies. We compare these phases to the conventional disordered Fermi liquid, and the 
insulating electron glass of semiconductors. Prospects for the realization of such glassy phases in cold atom 
systems are discussed. 
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I. INTRODUCTION 

There is much current interest in experiments with ultracold 
atoms and photons that provide clean realizations of models 
from condensed matter physics. A variant of the antiferro- 
magnetic Ising model in one dimension, for example, has re- 
cently been "quantum-simulated" with bosons in optical lat- 
tices 1 1 1 -an encouraging step toward quantum simulation of 
more general, strongly correlated quantum magnets in low di- 
mensions. The hope is that these quantum optics experiments 
eventually reach the parameter regimes and accuracy neces- 
sary to allow for predictions that can overcome the limitations 
of conventional theoretical approaches for strongly interacting 
quantum many-body systems. 

On top of that, the tunability of ultracold atoms allows one 
to explore new quantum many-body Hilbert spaces that have 
no direct condensed matter analog. In a series of remark- 
able experiments at ETH Zurich ||2}S|, Baumann et ai, have 
begun the quantum-simulation of strongly-interacting quan- 
tum gases with genuine long-range interactions [!5 |. In these 
many-body cavity QED systems, an atomic ensemble (a ther- 
mal cloud [6 7 ] or Bose-Einstein condensate) is loaded into an 
optical cavity containing quantized photon modes. Because 
the photons are massless, they mediate an interatomic interac- 
tion which does not decay as a function of distance between 
atoms and therefore couples all the particles in the ensemble 
to one another. Exploiting this property, Baumann et al. ||4]|8l 
found a mapping of their entangled atom-light system to the 
classic Dicke model describing two-level atoms uniformly 
coupled to a single quantized photon mode Ig UTTl . The su- 
perradiance transition of the Z2 Dicke model spontaneously 
breaks an Ising symmetry and may be viewed as a realization 
of an Ising ferromagnet which is exactly solvable due to the 
infinite range of the photon-mediated "spin-spin" interaction. 

Strack and Sachdev |12| recently computed the phase di- 
agram and spectral properties for atomic "spins" in multi- 
mode cavities assuming that the cavity mode functions and 
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(fixed) positions of the atoms can be chosen so that the effec- 
tive "spin-spin" interactions become random frustrated. It 
was shown that a quantum phase transition in the universality 
class of the (solvable) infinite-range Ising quantum spin glass 
lfT3lfT4ll occurs, potentially enabling comparison between ex- 
periment and the theory of spin glasses. In a related paper 
[15 -17], Gopalakrishnan et al. provided experimental details 
and detection methods for the spin glass phase. 

In the present paper we explore quantum glassiness in 
the charge or density sector of itinerant fermionic atoms in 
multi-mode cavities. An important difference to the previous 
study [12,1 consists in the inclusion of hopping of atoms on 
the lattice; the resulting phase diagram now depends on the 
quantum statistics of the itinerant particles. Related bosonic 
versions (Bose Hubbard models coupled to cavity photons) 
were considered in Ref. [191 wherein a superradiant Mott in- 
sulating phase, displaying entanglement of the charge of the 
atoms with a cavity mode, was found. Optomechanical ap- 
plications of degenerate fermions in a cavity were considered 
previously in Ref. |20] 

Our calculations provide evidence that multi-mode cavities 
with degenerate fermionic atoms can quantum-simulate vari- 
ous phases and properties of infinite-range glasses that share 
several properties with Efros-Shklovskii Coulomb glasses in 
the quantum regime. 

For a single-mode cavity coupled to itinerant bosonic atoms 
the onset of superradiance is concomitant with transla- 
tional symmetry breaking and the formation of a charge den- 
sity wave with a period of half the cavity photon wavelength. 
The glassy ordering of fermionic atoms in multi-mode cavities 
can instead be understood as the formation of one out of many 
energetically low-lying amorphous charge density patterns in 
a random linear combination of cavity modes. Most of these 
amorphous density orderings are metastable and the ensuing 
slow relaxation dynamics should in principle be measurable 
as a response to probe laser fields. 

Glassy phases of fermions have long been predicted to exist 
in Coulomb frustrated semiconductors, so-called electron- or 
Coulomb glasses fT\\\ in our case, the fermionic atoms play 
the role of the electrons in these earlier studies. In these elec- 
tronic systems frustration results naturally from a competition 
between long range Coulomb repulsions and the random po- 
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sitions and energies of impurity sites. While the former favor 
a regular density pattern, such as in a Wigner crystal, the latter 
disrupt this order and provoke a non-crystalline density order. 
The competition of these two ingredients leads to many long- 
lived metastable configurations with very slow relaxation dy- 
namics between them. This structure of phase space entails 
remarkable out-of-equilibrium properties such as memory ef- 
fects and logarithmic relaxation persisting over many hours 
|l22 -24| which are also present in quantum electron glasses 
approaching a metal insulator transition 1 25 i . 

An important hallmark of such insulating electron glasses 
is the Efros-Shklovskii gap in the single particle density of 
states in the glassy phase, whose elementary charge excita- 
tions are strongly Anderson localized. Such a pseudogap is 
required by the stability of metastable states in the presence 
of unscreened long range interactions [26 J . Remarkably, the 
amorphous ordering of the glass softens the hard gap -that 
would exist in a regularly ordered system- just to the maxi- 
mal extent which is still compatible with stability, leaving the 
glass in an interesting state of criticality Il27l|28l . This critical- 
ity results in a widely distributed response to a local excitation 
and avalanches |30|, features which occur also in mean field 
spin glasses |31 , 32 1. This criticality survives in the presence 
of weak quantum fluctuations (non-zero tunneling amplitude 
between sites in electron glasses or spontaneous spin-flips in- 
duced by a transverse field in spin glasses). It entails gapless 
collective excitations despite the absence of a broken contin- 
uous symmetry. Eventually the glass order melts at a critical 
value of the hopping or transverse field lfT4ll45ll50ll5TI . 

Upon further increasing the quantum fluctuations, a com- 
mon assumption is that the delocalization of the fermionic 
quasiparticles, i.e. the insulator-to-metal transition, co-incides 
with the disappearance of glassy dynamics, giving way to a 
disordered Fermi liquid. However, it is also possible that these 
two transitions are separated. The Anderson delocalization of 
the fermionic quasiparticles may precede the melting of the 
glass, in which case we obtain a metallic glass with non-zero 
conductivity at zero temperature. Such a glassy state with 
metallic conduction was obtained in dynamical mean field 
theory by Dobrosavljevic and collaborators 11331 [34l . More- 
over, the fact that glassy phases may also exist in phases with 
good transport properties was recently shown in models of 
frustrated bosons lfT8ll42ll43l where a 'superglass' phase with 
microscopically coexisting superfluid and glassy density or- 
der exists. In the context of Coulomb frustrated systems in 
condensed matter (without disorder), an intermediate metallic 
phase with periodic, striped density order ("conducting crys- 
tal") was discussed by Spivak and Kivelson l35l . 



A. Overview of key results and outline of paper 

In the present paper, we argue that two types of glassy 
phases, a metallic glass and an insulating Anderson-Efros- 
Shklovskii glass also exist for fermions with random infinite- 
range interactions (see Fig.[T]i. For low densities, we find that 
the metallic glass is avoided; instead, the liquid abruptly tran- 
sitions to the localized glass state (see Fig.|2]i. 



Our metallic glass state should not be confused with the 
metallic glass of metallurgy. In the latter materials, the glassy 
physics is entirely due to classical atoms freezing into out- 
of-equilibrium configurations, and the metalllic conduction is 
due to conduction electrons which move in the background of 
the frozen atoms. In contrast, in our system, the glassy dy- 
namics and metallic conduction are due to the same fermionic 
degrees of freedom, which are electrons in the condensed mat- 
ter realizations, and fermionic atoms in quantum optics real- 
izations. 

In Sec. [nj starting from a Jaynes-Cummings type Hamil- 
tonian for itinerant fermions coupled to cavity photons, we 
derive the fermionic model that we study in this paper: 

N ^ N 

(/.;> i=l /j=l 

(1) 

which contains a short-range hopping term t, disordered, 
random onsite energies e,, and long-range, random density- 
density interactions V/y mediated by photons. 

In Fig.[T] we show the phase diagram of this model at mod- 
erate density, n - (9(1), as a function of effective onsite dis- 
order (W defined in Eq. (26 1) and photon-mediated interac- 
tion strength (7) in units of hopping t. For small effective 
onsite disorder W > 0, as the interaction is increased, the 
disordered Fermi liquid (FL) becomes unstable to the forma- 
tion of an irregular, glassy density pattern, which depends on 
the interactions mediated by the random cavity modes. How- 
ever, the irregular density waves do not gap the Fermi surface, 
but leave the fermions metallic, with a finite conductivity in 
the low temperature limit. The glassy charge density order 
is marginally stable, which leads to soft collective density 
excitations. The scattering of fermions from collective den- 
sity modes leads to some non-Fermi liquid properties (such as 
finite-temperature transport) but the fermionic quasiparticles 
remain well-defined. The quantum glass transition and prop- 
erties of the metallic glass are analyzed with effective field 
theory methods in SecfV] 

Upon further reduction of the hopping, or increase of inter- 
actions, the random Hartree potential generated by the frozen 
density pattern starts to localize the quasiparticles and induces 
an Anderson insulator The latter has strongly suppressed dif- 
fusion, which is expected to vanish at T = 0. At this point, 
the metallic background vanishes, leaving behind an insulat- 
ing charge glass with a spatially strongly fluctuating frozen- 
in density distribution. This phase and an estimate for the 
transition point (//Oc.ioc for a 3t/ cubic lattice is presented in 
Sec. [Ill] We will show there that for low fermion densities, 
at equilibrium, the metallic glass is avoided and the Fermi 
liquid transitions discontinuously to a localized glass phase. 
Nevertheless, even at lower densities, the metallic glass may 
be experimentally observable, as it is expected to exist as a 
long-lived metastable phase, which eventually will nucleate 
the energetically more favorable localized glass. 

In Fig.[T] the metallic glass "strip" is predicted to separate 
the AES glass from the Fermi liquid even close to the Ander- 
son transition at ( J - 0, Wioc)- In that regime, as we explain in 
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FIG. 1: (Color online) Schematic phase diagram of model ^ for 
moderate densities in the vicinity of half-filling n a; 1/2 . y- is the 
variance of the long-range disorder V,^, and the variance of the ef- 
fective onsite disorder, Eq. \26) . The black, dashed line is the metal- 
insulator transition. The blue line is the glass transition. Quantitative 
computations in this paper are restricted to the W = Q axis, while the 
nature of the transition lines around the Anderson transition (J = 0) 
are inferred from scaling considerations in Sec. |VI| The localization 
transition from the metallic glass to the Anderson-Efros-Shklovskii 
glass and the signatures of this insulating glass phase are discussed 
in Sec.^l] We describe the Fermi liquid to metallic glass transition 
and the properties of the metallic glass phase in Sec.[V| 
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FIG. 2: (Color online) Schematic phase diagram of Eq. ([TJ for small 
effective onsite disorder (IV = 0) as a function of fermion density n. 
Below the crossing point (yellow dot) of the metal-insulator transi- 
tion (black-dashed line) and the glass transition (blue), the transition 
becomes first order. The Fermi liquid and the metallic glass may still 
exist (as long-lived metastable states) to the right of the dashed line. 



Sec. VI the quantum glass transition is sensitive to the fractal 
nature of the fermionic wavefunctions close to the Anderson 
transition. 



In Sec. VII we conclude with a unifying discussion about 
the potential of many-body cavity QED as a quantum sim- 
ulator of long-range quantum glasses exhibiting freezing in 
different sectors (spin, fermionic charge, bosonic charge). We 
also outline interesting open questions for future research. 



II. MODEL 

We consider spinless lattice fermions coupled to multiple 
cavity photon modes. The absorption of cavity photons (rep- 
resented by canonical bosonic creation and annihilation op- 
erators fl' , fl) raises the internal state of the fermions from 
a ground state (represented by canonical fermion creation and 
annihilation operators c^, c^) to an excited state (cj, Cg). In ad- 
dition, a classically-treated pump laser coherently drives tran- 
sitions between the ground state and the excited state. The 
Hamiltonian operator, 



^[Cgj Cg, a, C^,,C^] — //atom + //photon + //pump> (2) 

consists of three terms. The first one describes itinerant 
fermions on a d-dimensional (optical) lattice with / - I, ...,N 
sites with short-range, nearest-neighbor ((/, j)) hopping t 
(taken to be the same for ground and excited state); 



(ij) 

N N 
+ (f'> -1^ + ^) + X (^''« " ^) "'-^ 



(3) 



Here fi is the chemical potential, e,g, 6,e are the single- 
paiticle energies of the ground and excited state, respectively, 
and rii^g - ^c, ^ are density operators for the atoms in the 
ground state, and analogously for the excited state. Follow- 
ing Maschler et al. |44|, we write our equations in a frame 
rotating with the frequency of the pump laser and A describes 
the atom-pump detuning which can generally be chosen large 
A » f, e,>. 

We further have M cavity photon modes with frequencies 
a)( and spatially non-uniform atom-photon couplings g,f that 
generally depend on the atom's site ; and the characteristics of 
the cavity photon mode £. Finally, we have a pump term with 
amplitude /i, that does not involve photon operators: 

M N M 

//photon = ^ (Oea\ae - ^ ^ gie (c,' .c,- + c^^c,- ^aj) , 



I photon = ^ t^fflf 
N 

i=l 



= 1 t=\ 



2 ".a 



(4) 



The excited state can be adiabatically eliminated. In the afore- 
mentioned regime of large atom-pump detuning, one can ig- 
nore the dependence of the effective ground state interactions 
on the dynamics and spatial propagation of the excited state 
||44|| . Therefore, the e, ^ do not appear in the parameters of the 
resulting effective Hamiltonian, which reads: 

H[cl, Cg, fl] = - f ^ [Ag^j^s + I'^-c) + Yj ^'-s " ^ + TT "''g 
(/,;■> -=1 ^ ^ 



(iJ) 

M N M ^ 

+ Yi '^ealai + YYj ^X^"'-" ("^ "f) 
' ■ /=i e=i 



N M 

1=1 {m=l 



(5) 
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We drop the subscript g from now on. We include the expec- 
tation value of the term in the last line as a contribution to the 
single-particle energies. 



Zgitgim , t \ , , i 
—^{a^a„) + e; + ^ 



(6) 



In a final step, we integrate out the photons from Eq. (j5]l and 
get the expression: 

A' ^ N 



(7) 



The long-range density-density interaction written in a path 
integral representation 



M 



'i(gj(hihj cjf 



(8) 



makes the dependence on the bosonic Matsubara frequency 
Q. explicit fTTl. The magnitude of Vij is proportional to the 
amplitude of the driving laser /i, and can therefore be tuned 
flexibly. The sign and spatial dependence of Vij is determined 
by the choice of mode profiles of the cavity modes and pump 
lasers as well as the orientation of the lattice within the cavity. 

In this paper, we are primarily interested in the case where 
the photon mode functions gugjc in Eq. ([8]l can be realized as 
randomly varying in sign and magnitude in each disorder re- 
alization and with M sufficiently large, we assume the y,j(Q) 
to be Gaussian-distributed with variance 

6Vij{QL)6V,f{Q.') = [s,,5jf+6j,6ij)V{^,a)IN. (9) 

The overline represents a disorder average, and 6Vij is the 
variation from the mean value. Such a mean value only shifts 
the chemical potential and can be dropped. Further, we as- 
sume in this paper that couplings between different sites are 
uncorrelated. 

We note that it should also be possible to generate ran- 
dom and frustrated interactions of longer range by using other 
means than random cavity modes. For example one might 
employ a second fermion species to generate RKKY-type in- 
teractions among the primary fermion species. As in metallic 
spin glasses, such interactions decay as a power-law with dis- 
tance and oscillate with periods of the Fermi wavelength of 
the second species, which induces frustration. 

In the calculations below, the results only depend on the 
variance, which respects time-translation invariance 



y(Q, -O) s y-(0) 



(10) 



To capture the main effects, as in Ref.[T2l we may assume the 
simplified form: 



(11) 



where is a prototypical photon frequency representative 
of the spectral range of photons that mediate the inter-atomic 



interaction and v the disorder strength. For most of the paper 
we will concentrate on frequencies O «: (Wq and work with 
the static limit of the couplings 



6V^j6V,j, = {6,,6jj,+Sj,5u')j'/N, 
Wo ' 



J = 7(0) 



(12) 



The photon contribution to the single-particle energies in 
Eq. (j6|l and the pump term generate random local potentials 
for the fermions, which may additionally be superposed by a 
random lattice potential. We summarize all these effects by 
assuming random (Gaussian-distributed) onsite energies with 
independently tunable variance 



(13) 



III. 



INSULATING ANDERSON-EFROS-SHKLOVSKII 
GLASS 



In the absence of hopping, f = 0, the Hamiltonian 
(|7]i reduces to the classical Sherrington-Kirkpatrick (SK) 
model |56il of localized "spins" (describing presence or ab- 
sence of a particle). This spins are subject to a random longi- 
tudinal field e,, and kept at fixed "magnetization" M = 2n - 1 
where n is the fermion density. The low temperature glass 
phase of this model is understood in great detail. As illus- 
trated by the red graph in Fig.[3j a typical configuration of this 
glass exhibits a linear soft gap in the distribution 
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1^1 < I (14) 



of local Hartree fields 



dH ^ 



(15) 



with coeflicient a a; 4 ■ 0.31 = 1 .24 ijST] EH and Gaussian 
decay for \(p\ » J. The brackets (...) stand for the thermody- 
namic average. As compared to the canonical SK model an 
extra factor of 4 arises because we consider Ising degrees of 
freedom with magnitude s~. = n, - 1/2 = ±1/2. Remarkably, 
the soft gap ( 14 1 at low fields is universal, that is, independent 



of the strength of the random fields and the average magneti- 
zation (density) 1291 . A similar soft gap, the Efros-Shklovskii 
Coulomb gap ||26| is also known to exist in electron glasses 
with Coulomb interactions 



A. 



Anderson-Efros-Shklovskii (AES) glass with quantum 
fluctuations/finite hopping (/ 0) 



Upon turning on quantum fluctuations via a finite hopping 
t + Q, the fermions hop within the disorder potential of the 
above discussed Hartree fields. We refer to the eigenener- 
gies of the resulting single particle problem as single parti- 
cle excitations. In the limit f — > the latter are completely 
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FIG. 3: (Color online) Sketch of the evolution of the distribution of 
single particle energies p(E) as the hopping t increases. In the clas- 
sical Efros-Shklovskii limit, t = (red solid curve), p(E) coincides 
with the distribution of Hartree fields, which exhibits a linear pseu- 
dogap. In the localized AES glass at finite hopping (blue, dot-dashed 
curve) a pseudogap with vanishing density of states at the Fermi level 
p(0) = persists. A finite density of states at £ = emerges when 
the quasiparticles delocalize and the system turns metallic. The dip in 
the density of states at low energy gradually weakens as t increases. 



localized and their energies coincide with the local Hartree 
fields. As long as the hopping stays below a critical value 
(f < floe), the fermionic atoms remain Anderson-localized. 
In that regime the single particle density of states in a given 
local minimum of the glass has to vanish at the Fermi level 
(down to energies of order 1 / VA^, at T - 0), otherwise the 
state would be unstable with respect to charge rearrangements. 
This follows from arguments analogous to those for quantum 
Coulomb glasses 1261 f38l - l40l . As a consequence, in a typi- 
cal AES glass state, the linear compressibility (the analogue 
of zero-field-cooled susceptibility in spin glasses) vanishes, 
even though the full thermodynamic (field-cooled) compress- 
ibility is finite. Despite the vanishing linear compressibility, 
there is no hard charge gap, as there are charge excitations at 
any finite energy. The qualitative evolution of the distribution 
of single particle energies with increasing hopping is sketched 
inFig.|3] 

We now explain the key differences of the AES glass to 
other previously discussed glass phases and describe its trans- 
port properties. 

(Quantum) electron glass: The glassy electrons in strongly 
doped semiconductors are localized mostly due to quenched 
potential disorder, whereas in the case of the AES glass the po- 
tential is mostly interaction-induced and thus self-generated. 
If W could be made to vanish, the localization would be en- 
tirely induced by the random interactions in a glassy frozen 
state. Both, the quantum electron glasses and the AES glass 
have linearly suppressed density of states at low energies. 

Anderson or Fermi glass: In these localized insulators, 
where quenched external disorder dominates, the quasiparti- 
cles are localized despite the presence of weak short range 
interactions. The bosonic analogues of such insulators are 



often referred to as Base glasses ll36l [37l . Both fermionic 
and bosonic version of these glasses onsite-disorder dom- 
inated and do not exhibit glassiness in the sense of hav- 
ing many metastable configurations separated by high multi- 
particle baniers, except if strong interactions are present in 
addition to the disorder potential. Another difference is that 
the strong interactions in the AES glass cause gapless collec- 
tive modes, which are absent in less strongly correlated Fermi 
and Bose glass insulators. 

Mott glass: In fermionic lattice systems at commensurate 
filling with strong nearest-neighbor interactions and disorder, 
this Mott glass arises as an intermediate phase which is in- 
compressible due to a hard charge gap, but has a finite AC 
conductivity cr(cj) at any positive frequency 1411 . The main 
contribution to cr(w) comes from local particle-hole excita- 
tions, which become gapless due to disorder. The AES glass 
is very different from this Mott glass, in that it has no hard 
charge gap, exists at any filling and does not rely on external 
potential disorder. 

Ising spin glass in transverse field: If the occupation of 
sites n, = 0, 1 is thought of as an Ising variable, the model (j?]) 
looks similar to a mean field Ising spin glass in a transverse 
field r. The difference with the model considered in this paper 
consists mostly in the way in which quantum fluctuations act 
on top of the classical SK glass. This entails however a crucial 
difference of the respective phases to expect for weak long 
range interactions, as we will explain below. 

The transverse field SK model is well known to undergo 
a quantum glass transition from a glass to a quantum disor- 
dered state at a critical value F ~ J. The spins undergo glassy 
freezing when their interaction strength is bigger than the in- 
verse polarizability of the softest two level systems. One can 
show that the thermodynamic glass transition persists in the 
presence of random longitudinal fields (the equivalent of the 
disorder potential e, in Eq. (j?])), even though the latter spoils 
the Ising symmetry of the system. The glass transition then 
has the character of an Akneida-Thouless transition, like for 
mean field spin glasses in an external homogeneous or ran- 
dom field. The only symmetry which is broken at this transi- 
tion is the replica symmetry, which is rigorously established 
to occur in models with infinite-range interactions, while it 
remains controversial in finite-dimensional, short range inter- 
acting systems. 

If the Ising states of a spin are thought of as two states of 
a localized fermion (two-level systems), the transverse field 
SK model describes interacting, hut fully localized fermions. 
It follows from the above discussion that such systems have a 
non-glassy insulating phase. However, for our model a regime 
of localized but non-glassy fermions cannot exist. The reason 
is that there is no mechanism which gaps single fermion exci- 
tations away from zero energy, in contrast to the finite tunnel- 
ing amplitude F of two level systems which does gap the local 
excitations. Therefore, in the model considered here, the gap- 
less low-energy fermionic states are always unstable to glass 
formation for any small infinite range interactions Vjj (a sim- 
ilar argument wasput forward in Ref. [34l l. This is why the 
AES glass in Fig.lllabove Wioc extends all the way to the axis 
7 = 0. 
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Transport properties: The main contribution to the finite 
temperature transport in the AES glass comes from inelastic 
scattering of the quasiparticles, with an inelastic rate which is 
expected to decrease as a power law of T . At low T, where this 
rate becomes smaller than the level spacing in the localization 
volume of the fermions, the charge transport is expected to 
proceed first by power law hopping, and at lowest T by vari- 
able range hopping, assisted by inelastic processes among the 
fermions. This type of transport is drastically less efficient 
than in the metallic glass of Section [V] In the AES glass, the 
resistivity that should diverge with vanishing T . At the same 
time thermal transport may still be rather good due to collec- 
tive density modes, which remain delocalized. 

It is tempting to speculate that, upon advancing the tech- 
niques of Ref. 60 to many-body cavity QED, the next gener- 
ation of experiments could measure also the transport proper- 
ties of neutral, glassy fermions. 



IV. DELOCALIZATION TRANSITION OUT OF^THE 
ANDERSON-EFROS-SHKLOVSKII GLASS AT W « 

As mentioned above, in the localized, insulating phase the 
Efros-Shklovskii stability argument assures a vanishing den- 
sity of single particle states at the Fermi level. However, the 
latter becomes finite when the quasiparticles delocalize at the 
insulator-to-metal transition, where diff'usive transport behav- 
ior sets in. This takes place when the hopping becomes com- 
parable to the typical potential difference between a site and 
its neighbors. A precise calculation of this transition would 
require the solution of the glass problem including quantum 
fluctuations, to obtain the distribution of effective onsite dis- 
order, and the subsequent analysis of a delocalization transi- 
tion of quasiparticles. However, an estimate for the critical 
value (7)1 , below which the system behaves metallic, can 
be obtained as follows. 



A. Moderate density: n ~ 0(1) 



Even if W of Eq. (26 1 is negligible, the soft gap in the 
Hartree potentials acts like a rather strong onsite disorder for 
the fermions. To discuss delocalization at the Fermi level we 
need to consider sites with Hartree potentials of the order of 
the hopping e < t. At such energies, the density of states in 
the soft gap can be roughly approximated by the constant den- 
sity of states of a box-distributed onsite potential of width Wj 
given by 



From this we may infer an estimate for the delocalization tran- 
sition out of the insulating glass phase at 



4.6. 



/ c,loc 



(18) 



In this rough estimate we have neglected quantum fluctuations 
of the density order in the localized phase which are expected 
to increase the value of (//Oc.ioc further. This is because they 
weaken the density inhomogeneity and thus the onsite disor- 
der generated by the interactions. 



B. Low density: n « 1 

At low fermion densities, n <K 1, the pseudogap is irrel- 
evant for the delocalization transition because the gap is re- 
stricted to very small energies, far below the hopping strength 



needed for delocalization. The frozen fields ( 15 1 have a typi- 
cal magnitude Wj ~ ^lnJ. If this is the dominant contribution 
to disorder (i.e., for weak external disorder, W < Wj), the de- 
locaUzation instability arises in the glass phase when f ~ Wj, 
i.e., when the interactions are reduced below the value 



(-//Oioc ~ n 



-1/2 



(19) 



At larger external disorder, delocalization happens when t ~ 
W, independently of the interaction strength J. However, the 



interaction strength ( 19 1 is parametrically smaller than the in- 
stability iJ/t)met ~ rT^ of the Fermi liquid toward the metal- 
lic glass, which we will derive below in Eq. (|44]). This implies 
that at low densities, beyond the yellow, tri-critical points of 
Fig. |2] the localized glass state jumps discontinuously into a 
non-glassy Fermi liquid via a first order transition. The lo- 
cation of that transition can be estimated by considering the 
competition of kinetic energy cost ~ t and potential energy 
gain ~ ^lnJ of transforming the Fermi liquid into a fully lo- 
calized glass state. This yields 



(7/Ois. ~ n-'l^. 



(20) 



which scales in the same way with density as the delocaliza- 
tion instability of the localized phase, but is physically distinct 
from it. However, since the interactions are essentially infinite 
ranged, the system cannot simply nucleate a localized glass 
droplet locally, but has to undergo the localization transition 
more or less at once in the whole system. This suggests that 
the Fermi liquid phase is a very long lived metastable state, 
even far beyond the equilibrium transition point {J/t)i^t. 



W, 



1 



P{(p X t) at 



(16) 



For a box distribution of onsite disorder, Anderson delocaliza- 
tion of quasiparticles on a cubic lattice in c/ = 3 dimensions is 
known to occur at a critical disorder strength f531 



(-) 

\ t I 



16.54. 



(17) 



V. METALLIC GLASS AT W « 

In this section, we develop an effective field theory ap- 
proach for the Fermi liquid to metallic glass transition (the 
blue line for W ~ in Fig.[l]i and compute key properties of 
the metallic glass. We can write the path integral pendant to 
Eq. (j?]) in terms of collective Hubbard-Stratonovich fields for 
charge fluctuations p,(r) and Lagrange multipliers a, which 
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enforce p,(t) - c,(t)c,(t), where c, c are Grassmann variables 
representing the fermion fields and T will denote temperature. 
With this field content, the partition function corresponding to 

the model of Eq. (7i, Z = J DaDpDcDc e'^'^<'+'^"'\ is deter- 
mined by the action 



So^ I dT V c,(t) {dr - c,(t) 
Jo 

Xi/r 
dT ^ f (c,(T-)Cy(T) + c/t)c,(t)) , 

(ij) 

•^int = -- V I I drVijij -T')pi{T)pjiT') 

2 Jo Jo 
+ ^'^ (^/'^)'^;('^) - Pj('^)) + P Z'^)*^) ■ (21) 



In order to streamline notation, we will use 2 signs for all 
lattice site summations and integrations over imaginary time. 

We now integrate out the fermions, performing a cumulant 
expansion in the terms oc a to obtain an effective action for the 
density fields Z = J DaDpe-^^"-P^ with 

iJ,T,T' j,T 

-J]ajiT){{njiT)}o-pj{T)) 

+ 1 «/W«>(^')<n/(T)n>(T'))o + ... , (22) 



where <...)o denotes the quantum and thermal average with re- 
spect to the non-disordered free fermion action Sq. We intro- 
duce the auxiliary Hubbard-Stratonovich density fields p,- as 
an intermediate step to make the physics more transparent and 
to make contact with our previous work [12]. As in Ref.[T2l 
we will later integrate them out again to obtain an action for 
the order parameter fields Q"''. In App. |Aj we present an al- 
ternative route to derive the final self-consistency equations 
( 38pO i, which off'ers insight into the nature of the approxi- 
mations made below. 

We proceed by integrating over a, without keeping explicit 
track of the higher-order terms in a denoted by the dots in 
Below we drop those corrections. 



Eq. (22 1. Below we drop those corrections. However, the 
alternative approach in App. |A]can in principle resum them, 
at the expense of replacing the bare density-density correla- 
tor Eq. ( [24| with the full proper "polarizability" |13|. By 
dropping the interaction corrections to the latter, we operate 
at a level equivalent to the approximation used by Miller and 
Huse for the closely related infinite-range quantum spin glass 
problem lfT3l . Those authors obtained a good estimate for 
the quantum glass transition point when compared with more 
elaborate studies |54, 55 1 . 

It is convenient to express the resulting action in terms of lo- 
cal fluctuations in the onsite density dpiij) - piij) - {ni(T)}() = 



P, (t) - «: 

•,j,T,T' 



with the bare density-density correlator denoted by 

nil.' = <n/(T)«/T'))o . 
The effective onsite disorder consists of two terms 
£/ = Bi - Y{nj(T')}oVij(T - r') . 



(24) 



(25) 



Focusing on low frequencies we neglect the retardation, and 
find the disorder variance from Eqs. ([9][T3|) as: 



dsi dsj = 6ij (W^ + 7^(0)(«i)2) s 6ijW^ . (26) 

The dots in ( |23| ) stand for interactions between more than two 
6pi. 

We proceed by employing standard replica methods f46l to 
average over the disorder configurations of the y,y(Q) and e,. 
The resulting inter-replica interaction term ~ 6p* proportional 
to the variance of Vij is non-local in both, imaginary time and 
position space. We decouple this terms with an inter-replica 
matrix-valued Hubbard-Stratonovich field, that depends on 
two frequencies, Qf(D.,D.') ^ 6p';(D.)6p'l(n'), which is bi- 
local in imaginary time, but local in position space. Finally, 
we write the n-times replicated, disorder-averaged partition 
function Z" - J DQDp e with the action 

X ^Pnr)[-Tl'%..6p''j{r')- ^ ^'^P?(T)5pf (r) 



i,J,r,T',a 

i,j,n,n',a,b 

1 

2 



',a.h 



y(Q, Q') 



Qf{n,Q.')Qfi-^l,-Q.'} 



Qfi-n,-^l')6p''j(Q.)6p']iQ') 



(27) 



In the end we will take the replica limit « — » to ex- 
tract quenched averages. The dots stand for further replica- 
diagonal interactions between several 5p,, which are generated 



by the higher order cumulants in Eq. ( 22 1 



Incorporating local disorder into the polarizablity, which 
replaces n|°'^^, in a complete solution will actually generate 
replica off'-diagonal terms in 11 that self-consistently depend 
on W and J^{0){n}^. This significantly complicates the anal- 
ysis. As announced above, we restrict here formally to the 
W - Q slice of Fig. [l] and defer the quantitative analysis of 
onsite-disorder effects to future work |72|. 



A. Saddle-point free energy (W = 0) 



Without any truncation in the cumulant expansion ( 22 1, the 
functional integration would yield results for the correlators 
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of 6p that automatically obey the Pauli principle. However, 
after a truncation the latter might be violated. In order to 
proceed with a saddle-point analysis, we introduce into the 
action a global Lagrange multiplier A" conjugate to density 
fluctuations, enforcing that correlation functions still obey the 
Pauli principle on average. Namely, A" is determined self- 
consistently so that the equality 



(28) 



1=1 



is satisfied. The lefthand side can be viewed as the variance 
of the binomial probability distribution of lattice sites being 
either occupied (p,- - 1) with probability n or empty (p,- - 0) 
with probability 1 - «. We recall that n is the average density. 
n - Q corresponds to an empty band, n - 1 /2 to half-filling, 
and at n = 1 the band is filled. 

This global constraint is analogous to the "spherical" ap- 
proximation in spin glasses. It can be shown that relaxing the 
(exact) local constraints to a global constraint does not affect 
the quantum critical behavior 1 14 47 -49|. 



Noticing that the action in Eq. ( 27 i is quadratic in the p- 
variable, we can integrate it out exactly. Note that the bare 
density correlator 11^'^^, ■ ■ only depends on the differences in 
time (|r - t'|) and space (|/ - j\). The resulting action has a 
prefactor A^, the number of atoms, and thus, a saddle-point 
evaluation becomes exact. The free energy 



■ lim lim — In Z" 



(29) 



becomes a functional of the order parameter field 

= ^ Z ef^^' = ]^ J]m(r)Sp1(T')} (30) 

and the Lagrange multipliers A". Assuming that on the sad- 
dle point the latter take a replica-symmetric value A, the func- 
tional takes the form !F = Um„^o ji'^n with 

a,b,n 

T 



speaking unstable towards the breaking of replica symmetry. 
The latter signals the emergence of many pure states in the 
glass phase, the breakdown of full thermalization, i.e., ergod- 
icity breaking and the associated out-of-equilibrium phenom- 
ena at long time scales. As we show in App. A] the glass 
instability conditions obtained below (Eq. |38|43 1 indeed sig- 
nal the instability of a replica symmetric saddle-point toward 
replica symmetry breaking. 

For the Q"'' fields, the following ansatz is natural: 



(32) 



Here, the Edwards-Anderson order parameter qea shows up 
both in diagonal and off'-diagonal entries of Q"^. This ansatz 
in terms of qea and the (site- and disorder-averaged) dynamic 
density response D(Q) is the most general one, respecting 
replica symmetry and time-translational invariance. 

A non-zero value of q^A signals a frozen-in density distri- 
bution of the atoms: 

qEA = lim 4 y {dpi(t)6pf(0)) ■ 

In the glass phase, the spatially non-uniform on-site densi- 
ties differ randomly from their average value (depending on 
the state into which the glass freezes) and retain that value 
for infinitely long times. As in any glass, these frozen den- 
sity patterns are expected to depend sensitively on the details 
of the quench protocol or the preparation history in general. 
Note however, that only in the case of vanishing effective dis- 
order, W = 0, qEA serves as an order parameter for the glass 
transition. In the more realistic case of finite onsite disorder, 

9^ 0, the system still exhibits a thermodynamic glass tran- 
sition, as do mean field spin glasses in random fields, but the 
only symmetry to be broken in that case is the replica symme- 
try, since density inhomogeneities already exist in the Fermi 
liquid phase (cf. App.[A|. 

The average dynamic density response, 

D{^)^^Y,{6pt{£l)5pt{-^)) , (33) 

characterizes the response of the fermions to local, time- 
periodic modulations of the density. In terms of the 
parametrization ( 32 1, the free energy ( 31 1 obtains as 



+ 2 ([-n«»(Q,q)]"' + 2^)5"^ T ^-Yfrn^m^ + \jHmm^A - - «)« 



- J\Q)Q'''\-Q., il) 



(31) 



n 



where tr^b denotes the trace over replica indices and the dots 
stand for the neglected terms mentioned above. We have im- 
posed time-translational invariance of 2"*: 

Q"\^M) ^ e''''(Q,Q')fe+n',o/r. 

As long as we are primarily interested in the critical behav- 
ior and the properties of a single typical state close to the glass 
transition, we may proceed with a replica symmetric calcula- 
tion, even though the replica-symmetric saddle point is strictly 



In 



jHO)qEA 



+ 2A- r(D.)D(Q.) 



-n(")(Q, q)]"' +2A- JH0)D(0) 



(34) 



This free energy depends on the microscopic parameters of 
the original fermionic theory via the bare density response 

<«,(r)«/T'))o = r2 r i^.-'"(-'H/q(.,-.,)n(«)(Q,q), 

(35) 



9 



where n*"' is the real part of the particle-hole bubble, i.e. the 
convolution of two bare fermion propagators: 



rf^'Caq) = Re 



^ r d''k 
-T) ^Go(w + k + q)G()(w, k) 



Re 



A f(h)-f(M 



r d''k 

J (2^1 



with f(x) - l/(exp[ji;/r] + 1) the Fermi function and the 
fermion lattice dispersion. For a 3d cubic lattice, for instance 



- -2f (cos kjc + cos ky + cos fe,^ - fi , 



(37) 



with t the nearest-neighbor hopping matrix element and fi the 
chemical potential which fixes the lattice filling = for 
half-filling, that is, n = 1 /2). 



From the free energy functional (34i we will extract the 



phase boundary between the Fermi liquid and the metallic 
quantum charge glass and the associated quantum-critical dy- 



namics of the density correlations. Minimization of Eq. (34i 
with respect to qea, ^ and D(Q.) for each Q, yields a set of 
coupled saddle-point equations. The derivative with respect 
to D(Q) for f2 > together with the derivative with respect to 
qEA requires the density response to obey the selfconsistency 
relation 



D(Q) 



-I 



1 



(In)'' [-n(")(Q,q)]"' -H 2/1 - y2(Q)Z)(Q) 



(38) 



The lefthand side, when written as a geometric series, can be 
seen to express the self-consistent resummation of all cactus 
diagrams in the interactions Vjj. Minimization of Eq. (34i 



with respect to A gives back the global constraint on density 
fluctuations; 



ryD(Q) + ^EA = «(!-«) 



n 



(39) 



Finally, the minimization with respect to D{Q) determines the 
Edwards-Anderson parameter self-consistently: 



qEA 



r d\ 

J (2^(W 



/'(O)^EA 



(40) 



([-n('')(0, q)] ^ +2A- y2(0)D(0)) 
For vanishing interactions, 7(f2) = 7(0) = 0, we have ^ea = 



and Eqs. (38|39 1 have the free fermion solution: 



z)^°^(Q) = J -0[-n(°)(Q,q)] 



(41) 



Indeed, for free fermions, the constraint T YiQ^D^^\Q) 
n{\ - «) is automatically fulfilled, and thus A{J = 0) = 0. 



Note that Eq. (40 1 is a variant of the general glass instability 
condition: 



N 



1 



(42) 



where is the full density-density correlator. In Appendix [A] 
we present an alternative route to derive Eqs. ( 38|40|42 i. 



B. Numerical results for metallic glass: phase diagram and 
density response 

We now compute the glass transition line (Fig. |4| and the 
associated density response (Fig.|5]l, assuming no effective on- 
site disorder W = on a 3d-cubic lattice. We simultaneously 
together with the criticality condition de- 
(where it permits a solution with infinites- 



(36) solve Eqs. (OF 



rived from Eq. (|40 
imal ^ea): 



39 



"J (2^aI 



^'(0) 



([-n(0)(0,q)]-' + 2.l,,gi,,, - 72(0)D(0)) 



(43) 



This criticality has an important consequence for the dynamic 
response. By writing D(Q) - D{Q) - 6Dq, and expanding 
Eq. ( [38] l around Q = 0, we find that the condition (|43| ) entails 
a more singular low frequency behavior SDq^ ~ VI^^L as com- 
pared to the behavior of the bare density response 6D^^^ ~ |Q|. 
This is illustrated by the explicit solution of D{Q.) on a fre- 
quency grid in Fig. |5] Eq. (39i with finite ^ea. and more 
general arguments presented below ensure that this criticality 
extends into the glass phase. We will present an approximate 
analytical calculation of D(Q) at the glass transition and in the 



metallic glass phase in the following subsection V C 



For the numerical solution of these equations we discretized 
D(Q) on a frequency grid with varying step size up to a 
hundred grid points exploiting D(Q) = D(-Q). We em- 
ployed a modified version of Powell's Hybrid method algo- 
rithm for multi-dimensional root-solving ll52ll . For the 3d- 
numerical momentum integrations of the right-hand-sides, 
we employed the VEGAS Monte Carlo algorithm ||52ll . To 
avoid 6-dimensional integrations at each step of the multi- 
dimensional root solver, the momentum-integrated particle- 
hole bubble was catalogued as 4-dimensional array and then 
"quadm-linearly" interpolated in the integrands for the 3- 
dimensional q-integration (see App. |B] for some excerpts of 

nW(Q,^„^y,^,)). 

The resulting phase boundary between Fermi liquid and 
metallic glass is plotted in Fig. |4]for various fermion densi- 
ties. Along the phase boundary the glass transition is of con- 
tinuous nature, i.e. with a continuous onset of the Edwards- 
Anderson parameter ^ea (in the case W = 0). Like in quan- 
tum spin glasses, but in contrast to structural glasses (super- 
cooled liquids), the dynamic freezing and the thermodynamic 
glass transition are thus expected to coincide. Accordingly 
one expects the replica symmetry breaking pattern within the 
glass phase to be continuous (full replica symmetry breaking) 
which ensures criticality and the existence of gapless collec- 
tive modes [50] . 



One can see from Eq. (43 i that as a function of density. 



the glass instability of the Fermi liquid occurs when J ~ 
ilf(n)n^'^)-\ Using that n"(n) ~ n/Epin) ~ n^-^'''/t, we in- 
fer that in dimensions d = 3 the glass instability of the Fermi 
liquid scales as 



(J/t) 



,^ „-3/2+2/rf _ -5/6 
met " — " 



(44) 
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class of infinite-range, metallic quantum Ising spin glasses. 
Those are characterized by the interaction- and disorder- 
induced freezing of Ising degrees of freedom in the presence 
of a metallic charge sector with gapless fermionic excitations, 
which damp the order parameter fluctuations. This was orig- 
inally put forward in the context of metallic spin glasses by 
Sachdev, Read, and Oppermann fST), as well as by Sengupta 
and Georges |47|. Later on this universality class was further 
analyzed in the form of a Landau theory for a mean field ver- 
sion of the electron-glass transition out of the Fermi liquid to 
by Dalidovich and Dobrosavljevic [45 J . 



C. Dynamic density response in the metallic glass 



FIG. 4: (Color online) Numerically computed phase boundary of the 
metallic glass for various densities and W = 0. The crosses corre- 
spond to data points computed from a simultaneous solution of the 
saddle-point equations l |38|39|43"^ ; they are connected as a guide to 
the eye. Note that the phase boundary to the insulating glass is not 
plotted (see Fig.|2]for an illustration of both phase boundaries). 
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FIG. 5: (Color online) Comparison of the collective density response 
as a function of imaginary frequencies at the glass transition (blue 
crosses on the frequency grid computed from saddle-point equa- 
tions l |38|39|43") , and its bare pendant from Eq. l |41| l. The slope of 
the low energy response in the glass phase is singular, with a Vl^ 
behavior for small frequencies. 



It is instructive to derive an approximate analytical solu- 
tion for D(Q) neglecting the q-dependence of the particle- 
hole bubble in Eq. To this end we take n<°>(Q, q) -» 
n*°\Q, Q) with a fixed, prototypical Q diff'erent from any 



potential nesting vectors, |Q| < 2|kp|. Then, Eq. (38 1 becomes 
a quadratic equation for the (approximate) density response, 
which we call Dq(Q.). It can be solved in closed form; 



£>q(Q) 



-n('')(Q,Q)]" 



+ 2A 



2JHQ) 



('[-n«')(Q,Q)] 



y2(Q) 



(45) 



where we chose the physical solution with minus-sign 
in front of the square-root, ensuring that Dq(Q.) decays 
to zero for large frequencies. Recall that the particle- 
hole bubble becomes zero for large external frequencies: 
liminHM n(°>(Q, Q) -> 0. The ciiticality condition, Eq. (|43|) 



for fixed momentum Q, applied to Eq. (45 1, yields: /lc,giass = 

7(0) - i [-n(0)(0,Q)]"', and DQ(0),,giass = jjoy At the glass 
transition Z)q(Q) develops a singular response to small fre- 
quency pertubations of the local fermion density. The approx- 
imate energy scale 



A + 



-n(0)(o,Q)] 



1\ 



(46) 



at low densities. This scaling with the fermion density is 
confirmed by the numerical results of Fig. |4] However, as 
we argued in Sec. IV B at lower densities (beyond the yel- 
low tricritical points in Fig. [T]), this instability is pre-empted 
in equilibrium by a first order transition which takes place 
at (7/f)ist ~ «"'^^ <& (//Omet- Nevertheless, the Fermi liq- 
uid phase should remain experimentally relevant even at these 
low densities, since it is very difficult to nucleate the localized 
insulator out of the metastable Fermi liquid phase, owing to 
the high nucleation barrier imposed by the long range interac- 
tions. 

The glass transition, Fig.|4] and the associated emergence of 
a singular density response, Fig.|5j are that of the universality 



controls the distance to the transition and separates two qual- 
itatively distinct regimes in the dynamic density response of 
the Fermi liquid. 



Continuing Eq. (45 i to real frequencies, iQ. ^ Q. + iO+, we 
obtain: 



-lmDQ(Q) 



-V^, Q» A, 
n, Q « A . 



(47) 



The self-consistency condition ( 39 1 pins A = in the en- 



tire glass phase, and the density response remains singular at 
zero frequency. This property holds true independently of the 
above approximations, as can be formally derived assuming 
full replica symmetry breaking in the glass phase 1/501 . 
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The key ingredient for the unusual square-root behavior for 
both, Dq(Q.) and the full D(Q) of Fig. [5] is the coupling of 
the collective charge fluctuations to the metallic particle-hole 
continuum which coexists with glassy, amorphous density or- 
der The low frequency behavior of the particle-hole bubble 
(see Fig.[6]in the appendix) at finite momentum transfer Q: 

rf°^(n,Q) -^"'(0,0) 



\n\ 



(48) 



underlies Eq.j47]i and also enters into the numerical compu- 
tation of Fig. p] 



D. Fermionic quasiparticles in the metallic glass 

The single-particle properties of the underlying fermions in 
systems belonging to the universality class of infinite-range, 
metallic quantum Ising glasses have been worked out by Sen- 
gupta and Georges BtI . At T - Q, at the critical point and 
in the glass phase, the fermions remain well-defined quasi- 
particles. Indeed, the leading self-energy correction due to 
low frequency charge fluctuations scale as 



3/2 



(49) 



with a frequency exponent > 1. Nevertheless, at finite tem- 
peratures this translates into non-analytic corrections in the 
transport properties. The finite-temperature resistivity, for ex- 
ample, scales as 6p(T) ~ T^^^ in the quantum-critical regime 
above the QCP |47| and in the entire metallic quantum glass 
phase. Deeper in the glass phase, the metallic diffusion even- 
tually breaks down when the localization transition to the 



Anderson-Efros-Shklovskii glass of Sec. Ill is reached 



VI. PHASE BOUNDARIES OF THE METALLIC GLASS AT 

W>J 

In the preceding section we have argued for a metallic glass 
phase at moderate densities and weak external disorder How- 
ever, also when the eff'ective onsite disorder W is larger than 
the infinite-range interaction J, we expect an intermediate 
metallic glass phase "strip" between the Fermi liquid and the 
AES glass, as shown in Fig.[T] Below we will present scaling 
arguments to justify this scenario, with both transitions at the 
borders of the metallic glass "strip" (the metal-insulator tran- 
sition on "top" and the glass transition on the "bottom") being 
continuous. 

For vanishing interactions 7 = 0, the two instability lines 
must join at the Anderson transition of the disordered free 
fermions, where Wioc/t - W\oc/t = 21.29 (for Gaussian dis- 
order at half filling f531). In the presence of weak interac- 
tions, the frozen fields from the glassy density order tend to 
increase the disorder variance by 6{W^) < J^, which has a 
similar effect as a weak increase of disorder 5W ~ J^/W. Ac- 
cordingly one expects the critical value W/t for delocalization 
to decrease by a relative amount 



w 




w 


(0) 


W 


if) 


_ t . 
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. t . 




. t . 


loc 



\i1 



(50) 



since stronger hopping is necessary to compensate for the ex- 
tra disorder. 



A. The role of fractality for the glass transition 

On the other hand, we have to analyze the glass instability 
( |42l ) within the metallic phase. In the limit 7 — » 0, it holds that 
X — > xi J - 0), which reduces to the (exact) susceptibility of 
non-interacting fermions. As noted in Ref. [34l the sum over 
susceptibility squai^es. 



- ^ - 0) 



(51) 



diverges at the Anderson transition. This implies a glass in- 
stability already within the metallic phase, even for very weak 
interactions 7 <K f, W. To the best of our knowledge, the pre- 
cise divergence ofxi^ or equivalently, o^x]p [t"] ~* [t"]! *' 
is not known. However, the disorder-averaged density-density 
correlator has been well studied, since it reveals interesting 
properties of the fractal nature of the electronic wavefunctions 
and the anomalous diffusion at the Anderson transition |62|. 
In particular, the spatial Fourier transform of Xiji^ - 0) be- 
haves as 



1 



(52) 



where ^ is the correlation length, which diverges at the Ander- 
son transition as ^ ~ (1 -f/fc)"^; 6^ - 1/v^"' is the single parti- 
cle level spacing in the correlation volume and v is the density 
of states. The scaling function fix) behaves as 116111631 



fix) 



X » 1, 

X <sc 1. 



(53) 
(54) 



Here, d2 is the fractal dimension associated with the distri- 
bution of |iA^'*(x)| over all space, iff(x) being a critical single 
particle wavefunction at the delocalization transition. From 
this one obtains that 



X2 



N ■ 
•J 

■2(d-d2) 



(0 = 0) 



(i-t/t,)- 



r d^q — 

J (2^0^^' 



0) 



(55) 



The exponent results as a = 2v(d - dz) ~ 4, using the known 
values d2[d^3]^13 JgH |6l and y = 1.57. [33l 

It is reasonable to expect that;^f2 ~ (1 - t/tc)^" needed in 
Eq. (51 1 diverges at least as fast SiSX2, and thus a > a. Hence 



we expect that the glass instability is displaced from the non- 
interacting Anderson transition by an amount 



(56) 
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Comparing with Eq. (50i and noting that a > 1, one sees 
that as 7 — » the glass instability line approaches the non- 
interacting Anderson transition from smaller values of (W/t) 
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than the delocalization instability, as indicated in Fig.[T] These 
arguments confirm the existence of an intermediate phase be- 
tween disordered Fermi liquid and Anderson localized glass, 
which exhibits both metallic diffusion and glassy density or- 
der. 



Vn. CONCLUSION 

Based on the calculations and arguments of this paper and 
a previous work I12J , we believe that many -body cavity QED 
could evolve into a new platform to explore the physics of 
long-range quantum glasses. The long range of the photon- 
mediated interactions simplifies the theoretical analysis for 
these systems; this may allow for quantitative comparison be- 
tween experiment and theory. Our work complements previ- 
ous proposals on glassy many-body systems in quantum op- 
tics ll69l iTOl which have focussed onto the creation of onsite 
disorder or random short range exchanges by using random 
optical lattices and species admixture. 

In Ref.[T2l we predicted a quantum spin glass transition of 
fixed, stationary atoms where the source of quantum fluctua- 
tions was spontaneous tunneling between suitable chosen in- 
ternal states of the atoms. In the present paper, we considered 
itinerant fermions where the source of quantum fluctuations is 
tunneling between adjacent lattice sites. We found a metallic 
glass phase with a gapless Fermi sea in the presence of ran- 
dom density order and, for stronger interactions, a localized 
state with strongly random charge distributions and vanishing 
conductivity at T = 0. 

Various aspects of the phase transitions into and out of these 
glass phases are worth further studies: One is to quantify fur- 
ther the impact of the fractality of wavefunctions close to the 
Anderson transition and compute the critical exponent a with 
which the glass instability line approaches the non-interacting 
Anderson transition point. Experimentally, repeated ramps 
from the Fermi liquid phase into the glass phase should pro- 
duce metastable states with a different density pattern at each 
run. The resulting variations in absorption images could be 
particularly strong, and therefore perhaps easier to detect, in 
the low density regime where the glass and localization tran- 
sition is of first order. 

What happens when one considers models of type Eq. ([T]i 
with bosonic atoms? Numerical simulations |42], '43 1 and 
replica calculations |18| for similar models agree on the ex- 
istence of a stable "superglass" phase with superfluid phase 
coherence in the presence of glassy random density order. It 
is tempting to identify the superglass as the bosonic pendant 
to the metallic fermion glass. However, at finite temperatures, 
the domain of the superglass (cf. Fig. 1 of Ref. 18) is reduced 
because the main effect of thermal fluctuations is to weaken 
the superfluid phase coherence of the bosons. For fermions, 
on the other hand, the domain of the metallic glass increases 
at finite T because the main effect of thermal fluctuations is to 
enhance the inelastic scattering rate and to weaken the local- 
izing disorder potential. 

We described the disordered atom-light quantum phases 
using an effective equilibrium ground-state description. In 



actual optical cavities, one deals with pumped, steady-state 
phases certain features of which may not be captured in an 
effective equilibrium description IfTTl |66l l67l . We hope to 
address the non-equilibrium properties of open, disordered 
Dicke models in the near future in a forthcoming paper. We 
also want to analyze how the intriguing properties of classical 
and quantum glasses such as aging (see Ref. 68 and references 
therein), or out of equilibrium dynamics and avalanches OTII 
are modified in the open, driven, steady-states in many-body 
cavity QED. 

Finally, it will be interesting to investigate how the results 
of this paper and Ref. 12 scale with the number of cavity 
modes and the choice of cavity mode profiles. 
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Appendix A: Alternative derivation of the glass instability 
equations j38|43|42> 

1. Replica approach 

Here we formally derive the glass instability for arbitrary 
onsite disorder using the replica approach. We start from 
model (j7|i. Replicating n times and taking the disorder av- 
erage and introducing Hubbard Stratonovich fields to decou- 
ple the quartic density interactions, one obtains the replicated 
partition function 

r Y\Dc'lir)Dc1(T) Yl DQab{T,T')exp[-Sa-S,] 

a,i,T a<b,T,T' 

(Al) 

with the action 

n „i/T N 
Sc^J] dTj]-c1iT)(dr-lu)c1iT) 

fl=l ^" i=I 

+ J f'^'^^Tj' (^"(^)^;^^) + c'jirXir)) 

a=\ (,-,;> 
1 " ni/T 

- 2 Z + -^'^^ " r')Q,t(T, t')] n1{T)n1{T') , 

a,b=l ^0 

= T Z [' dTdT'j\T - t')QI,(t, t') . (A2) 

In the thermodynamic limit (A^ oo), the infinite range of 
the interactions allows us to take the saddle point with respect 
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to Qab, which satisfy the saddle point equations 

ai(T,T')= ^^K(T)«f(T')). 



(A3) 



(Note the slightly different definition of Q as compared to Qa 



defined in Eq. (30i.) From here on we neglect the retardation 
in the mediated coupling J(t - t') and replace it simply with 
J. As usual, the saddle point values of Qai, are independent 
of T, t' for a b, and depend only on r - t' for the replica 
diagonal Qa^. In the disordered, replica symmetric phase the 
term = + J^Qai^t can be recognized as the variance 
of the self-consistent effective disorder p5| ). Note that since 
Qai^b ^ (n)^, the effective disorder never really vanishes, un- 
less e, and y,y have special correlations. 

To detect a glass instability one has to solve first the saddle 
point equations of the disordered phase. This yields a replica 
symmetric solution G^fir, t') which encodes the Edwards An- 
derson overlap 



(A4) 



and the average local susceptibility, 

<3aa ("T - ■r') - ?EA = ^<n;(T)«/(T'))c = ^1oc(t - t') . 

(A5) 

The glass instability is found by writing <3 = GF^ + 6Q 
and expanding the free energy in 6Q. A standard cumulant 
expansion yields 

Z^[Q^^ + 6Q] = exp[-N(fiF^^ + 6(fiF))] , (A6) 

which by virtue of the extremality of Q^^ starts with a 
quadratic term. The glass instability is signalled by the vanish- 
ing of the coefficient of the term <52^^ (the mass of the replicon 
mode). 



= (/3JY 



N 



dTdT'niJjjrij^aiT') 



1 - 



'J 



4 



(A7) 



Hereby the correlator 

Xij^ j dT{nUr)njAO))^ , (A8) 

has to be evaluated with the replica symmetric action. 
Eq. ( A7 1 is to be compared with Eqs. ( 42|43 i. 



2. Cavity approach: local selfconsistent action 

In this subsection we derive the glass instability conditions 
from a selfconsistent local action derived within a cavity ap- 



particle and interaction part S = Si + Sy, 



pin 

Sx^ \ dT) Ciir) (dr + Si - ju) c,(t) 

Jo tr 

dr t (c,(t)c/t) + c/t)c,(t)) 



-X 



dT'Vn{T-T')X 



(«,(t) - <n/»(n/T') - {nj)) , (A9) 
where ni{T) = Ci{T)ci{T), and 

Si = Si-J^V,j(n^O){nj). (AlO) 



We now take advantage of the infinite range nature of the in- 
teractions, to transform the above problem exactly into a self- 
consistent single-site problem with retarded density-density 
interactions. The extra contribution to the disorder ( |A10| l can 
be treated as an additional Gaussian disorder with variance 
This type of disorder is essentially 

unavoidable in the system. Note that if it does not vanish, 
e, 0, there are density inhomogeneities already in the dis- 
ordered phase, i.e., («,) {rij) for / j. Unfortunately this 
renders the exact evaluation of the self-consistency problem 
very hard, and one has to resort to approximations in order to 
obtain quantitative results. 

Upon average over Vij a subsequent Hubbard Stratonovich 
transformation and a saddle point approximation, one finds 
that the interaction term can be resummed as 

S'y = ^(«,(t) - {ni))R{T - t')(«,(t') - <«,», (All) 
with the kernel 



{All) 



where the average local susceptibility ;\fioc must be found self- 
consistently. 

When computing the density-density correlator 



;e,/Q) = <n,(n)n>(-f2)>5,+5',, 



(A13) 



with the self-consistent action, one should be aware, however, 
that the above resummation does not include terms which con- 
tain a given coupling Vij only once within the expansion in 
interactions. Those are indeed irrelevant upon site or disorder 
averaging. However, they cannot be neglected when the glass 
susceptibility is computed. 



proach. We start from the action (21 1, and split it into single 



(A14) 
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A glass transition manifests itself by a divergence of this sus- 
ceptibility, which in mean field systems signals the emergence 
of ergodicity breaking in the form of many pure states, and 
replica symmetry breaking. In this formula Xij is the full 
density-density correlator. To the relevant order in Vij it can be 
obtained from Xjj by the summation of the geometric series. 



-■Xij + 



(A15) 



Upon taking the square and summing over all pairs of sites 
one finds 



1 



Xglass — 



(A16) 



which diverges when the denominator vanishes, reproducing 



Eq. ( |A7 i. 

In mean field glasses, the continuous breaking of replica 



symmetry ensures usually that the condition (A16l remains 



fulfilled, even within the glass phase, where Xij is to be inter- 
preted as the density-density correlator within one metastable 
state (and contributions from single couplings dropped). This 
was explicitly shown in the case of the infinite range trans- 
verse field Ising spin glass |50|. This phenomenon of main- 
tained marginal stability is at the basis of the permanent gap- 
lessness of the quantum glass phase. 



3. Generalized Miller-Huse type analysis 

Following a reasoning similar to the one by Miller and 
Huse [TTl for the transverse field spin glass, we formally com- 
pute the local density-density correlator in a perturbation se- 
ries in the interaction action Sy- The perturbation series for 
the correlator ;f- can be formally summed up as a geometric se- 
ries of local interactions R(t - t') linking "proper polarizabil- 
ity" blobs n,j(T - t') that cannot be reduced into two separate 
pieces by cutting a single /?-line. 



i„(Q) = n,/Q) + 2 ntt(Q)/?(Q)n,/Q) + 

k 

1 



(A17) 
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The proper polarizability has itself a power series expansion 
in R. To lowest order one has simply 



n,/Q)=;r(°'(Q) + OW, 



(A18) 



where denotes the non-interacting density-density corre- 
lator. 



(A19) 



As derived above in Eqs. ( A7|A16 1, the glass transition oc- 
curs when 



1 ^ :^^Tr [n '(f^ = 0) - 7?(Q = 0)]" 



N 



(A20) 



In the approximation where we neglect the effective disorder, 
e, = 0, in the quantum disordered phase, one has translational 
invariance, which allows one to work in Fourier space. 
The average local density-density correlator, 

^(O) = ^ Yj>^im = ^Tr [n-'(n) - , (A21) 



must satisfy the selfconsistency condition 
R{D) = 7^(Q)D(n). 



(A22) 



It must also obey the constraint 

J ^Z)(Q)= i^(«,(T)-<«,»2 (A23) 



In the absence of onsite disorder, the last term in brackets 
vanishes. The constraint (A23i is fulfilled automatically by 
an exact solution. 

" one 



However, if D is evaluated within an ap- 
proximate scheme, e.g. with the help of Eqs. ( A17|A18 i 



should impose this constraint to obtain a better approximation. 



In particular, we can satisfy the short time constraint ( A23 i by 
adding an adjustable short-time component A to the relation 
J^(Q.){D(Q) - A) - R(D.), to correct for the errors at high fre- 
quencies introduced by the approximations involved in com- 
puting D. The better the approximation, the smaller will be 
the A required to enforce the short time constraint. This recipe 
turns out to be essentially equivalent to the global constraint 
we introduced in Sec. IV A| 

We thus have to solve simultaneously 

J ^D{n)^n{l-n)-[(^-n\ (A24) 



and 



Z)(Q) 



T I 



1 



y2(Q)[D(Q) - A] 



(A25) 



The glass transition arises when 



^(0)^ 



Tr 



1 



\n-i(0)-y2(0)[D(0)-^] 



d 



1 



-Tr 



1 



d[D(0)]N \n-HO) - JH0)[D(0) - A] 



(A26) 



The latter relation leads to a singular behavior of D(Q.) around 
Q — » 0, as one may see by expanding Eq. ( |A25| l around Q = 
0, very similarly as in quantum spin glasses [13, 50 1 . This 
singularity ensures the presence of spectral weight ImD(Q — > 
w + id) at all finite real frequencies oj. 



The three last equations are to be compared with Eqs. (38 



40 1 to which they reduce in the translationally invariant case. 
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Appendix B: Particle-hole bubble as function of external 
frequency and momenta 

In Fig. |6] we display the particle-hole bubble as a function 
of external frequency for fixed momentum transfer. The low- 



frequency part behaves as \Q.\ as alluded to in Eq. (48 i. 

In Fig. |7] we plot an exemplary Q = bubble as func- 
tion of momenta as occurring in the saddle-point equations 
( 38|39|43| ). As expected, the dominant contributions come 
from momenta in the vicinity of the nesting condition q ^ 
Quest = i^,^,^)- Although logarithmically divergent at q = 
Qnest, the right-hand-side of the saddle-point equations re- 
mains regular as it involves an additional 3 -dimensional in- 
tegration over q. 



FIG. 6: (Color online) \Q.\ behavior of the particle-hole bubble, 
Eq. \36\ , as a function of external frequency for fixed momentum 
transfer Q. 




FIG. 7: (Color online) Exemplary momentum behavior of the static 
(Q. = 0) particle-hole bubble, Eq. l|36|l, at half-filling in 3 dimensions. 
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